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The present surge for the astrophysical relevance of boson stars stems from the speculative 
possibility that these compact objects could provide a considerable fraction of the non-baryonic 
part of dark matter within the halo of galaxies. For a very light 'universal' axion of effective string 
models, their total gravitational mass will be in the most likely range of ~ 0.5 M of MACHOs. 
According to this framework, gravitational microlensing is indirectly "weighing" the axion mass, 
resulting in ~ 10~ 10 eV/c 2 . This conclusion is not changing much, if we use a dilaton type self- 
interaction for the bosons. Moreover, we review their formation, rotation and stability as likely 
candidates of astrophysical importance. 
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I. INTRODUCTION 
A. Dark matter — Issue of missing mass 

The rotation velocities of spiral galaxies can be accurately measured from the Doppler effect. At large radii where 
the stellar surface brightness is falling exponentially, velocities are obtained for clouds of neutral hydrogen using the 
21 cm hyperfine line. The resulting 'rotation curves' are found to be roughly flat out to the maximum observed 
radii of about 50 kpc, which implies an enclosed mass increasing linearly with radius. This mass profile is much 
more extended than the distribution of starlight, which typically converges within ~ 10 kpc; thus, the galaxies are 
presumed to be surrounded by extended halos of dark matter. 

Perhaps the most compelling evidence for dark matter comes from clusters of galaxies. These are structures of 
about 1 Mpc size containing more than 100 galaxies, representing an overdensity of about a factor 1000 relative to 
the mean galaxy density. It is assumed that they are gravitationally bound since the time for galaxies to cross the 
cluster lasts only about 10% of the age of the Universe. The cluster masses are determined in several independent 
ways: First, the virial theorem uses the radial velocities of individual galaxies as 'test particles'. Second, observations 
of hot gas at about 10 7 K contained in the clusters, which is observed in X-rays via thermal bremsstrahlung. The gas 
temperature is derived from the X-ray spectrum, and the density profile from the map of the X-ray surface brightness. 
Assuming the gas is pressure-supported against the gravitational potential leads to a mass profile for the cluster. 
The third method is gravitational lcnsing of background objects by the cluster potential. There are two regimes: the 
'strong lensing' regime at small radii, which leads to arcs and multiple images, and the weak lensing regime at large 
radii, which causes background galaxies to be preferentially stretched in the tangential direction. All three methods 
yield estimates for cluster masses pL^ which show that visible stars contribute only a few percent of the observed 
mass, and the hot X-ray gas only about 10-20%, hence, clusters are dominated by dark matter. 

On the largest scales, there is further evidence for dark matter: 'streaming motions' of galaxies (e.g., towards nearby 
superclusters such as the "Great Attractor" ) can be compared to maps of the galaxy density from redshift surveys 
to yield estimates of f2. Here the theory is more straightforward since the density perturbations are still in the linear 
regime, but the observations are less secure. A similar estimate may be derived by comparing our Galaxy's 600 km/s 
motion towards the Virgo cluster relative to the cosmic rest frame, confirmed by the observed temperature dipole in 
the cosmic microwave background (CMB). 
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Furthermore, it is possible to connect the observed large-scale structure in the galaxy distribution with the results of 
the CMB anisotropics if the universe is dominated by non-baryonic dark matter. Commonly, the present matter/energy 
density £1 = Qm + S1a of the Universe is decomposed into two components. There is accumulating evidence for 
ilo = 1 ±0.2 and (total) matter density S1m = 0.4 ±0.1 which implies a vacuum energy or dimensionless cosmological 
constant of f^A = 0.85 ± 0.2 cf. M- Theories of inflation prefer a flat Universe with f2 = 1 as its most 'natural' 
value; this also requires non-baryonic dark matter. 

B. Dark matter — Candidates 

What are realistic candidates for dark matter? Hot gas appears to be excluded Q| by limits on the Compton 
distortion of the blackbody CMB spectrum; atomic hydrogen due to 21 cm observations; and ordinary stars by faint 
star counts. Asteroids are very unlikely since stars do not process hydrogen into heavy elements very efficiently and 
hydrogen 'snowballs' should evaporate or lead to excessive cratering on the Moon. Black holes more massive than 
~ 10 5 Mq would destroy small globular clusters by tidal effects. 

Today's most viable dark matter candidates fall into two broad classes: astrophysical size objects called Missive 
Compact Halo Objects (MACHOs), and so-called Weakly Interacting Massive Particles (WIMPs). Actually these 
classes could possibly be interrelated, as we are going to show. 

Several different objects belong to the first class: Jupiter-size brown dwarfs consisting of hydrogen and helium less 
massive than 0.08 Mq are the most prominent possibility. Below this limit, the central temperature is not sufficient 
in order to ignite hydrogen fusion, so these objects just radiate very weakly in the infrared due to gravitational 
contraction. Other MACHO candidates include stellar remnants such as cool white dwarfs, neutron stars, and 
primordial black holes or as a result of a supernova. 

The WIMP candidates are the invisible axion (hypothesized to solve the strong CP problem or reemerging 'univer- 
sally' in effective string Lagrangians), one of the neutrinos (provided it has a mass of about 10 eV), and the lightest 
supersymmetric particle, the neutralino, which is expected to be stable. All these have to have a very weak interaction, 
so that they could not be detected so far. 

C. Gravitational microlensing 

The conclusion of gravitational microlensing of stars within the Large Magellanic Cloud (LMC) is that MACHOs 
in the planetary mass range 10 -6 to 0.05 M© do not contribute a substantial fraction of the Galactic dark halo. In 
the two-year data of the LMC events of the MACHO group 8 events could be detected which are well in excess of 
the predicted background of approximately 1.1 events arising from known stellar populations. Hence, MACHOs in a 
dark halo appear to be a natural explanation. 

A statistical analysis of the galactic halo via microlensing HJ| suggests that MACHOs account for a significant 
part (> 20%) of the total halo mass of our galaxy. Their most likely mass range seems to be in the range 0.3 - 0.8 
Mq, with an average mass of 0.5 M , cf. g(Jg. If the bulge is more massive than the standard halo model assumes, 
the average MACHO mass |l(| will be somewhat lower at ~ 0.1 Mq. 

This can be viewed as an indication that MACHOs form an distinct large class of old objects that cannot be easily 
extrapolated from any familiar stellar population, such as brown or white dwarfs. 

However, there are some astrophysical difficulties with this interpretation, mainly arising from the estimated mass 
~ 0.5 Mq for the lenses. These cannot be hydrogen-burning stars in the halo since such objects are limited to less than 
3% of the halo mass by deep star counts (TTJ . Modifying the halo model to slow down the lens velocities can reduce 
the implied lens mass somewhat, but probably not below the substellar limit 0.08 M . Old white dwarfs have about 
the right mass and can evade the direct-detection constraints, but it is difficult to form them with high efficiency, and 
there may be problems with overproduction of metals and overproduction of light at high redshifts from the luminous 
stars which were the progenitors of the white dwarfs ]l2|| . Primordial black holes are a viable possibility, though one 
has to appeal to a coincidence to have them in a stellar mass range. 

Due to these difficulties of getting MACHOs in the inferred mass range without violating other constraints, there 
have been a number of suggestions for explaining the LMC events without recourse to a dark population: most of 
these suggestions construct some non-standard distribution of 'ordinary' stars along the LMC line of sight. However, 
these proposals appear somewhat contrived 0], but can be tested observationally in the near future. 
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D. Boson stars or axion stars as alternatives? 

It has been recently suggested jl3| that MACHOs could rather be primordial black holes formed during the early 
QCD epoch in the inflationary scenario. For cosmological dark matter, bound states of gravitational waves, so-called 
'gravitational geons' built from spin-2 bosons, are also considered recently Q. 

Since the standard model of elementary particles as well as their superstring extensions involve also Higgs type 
scalar fields, we will analyze here the alternative possibility that primordial boson stars account for this non-baryonic 
part of dark matter Jl5||l6| . Boson stars are descendants of the so-called geons of Wheeler [|l7|— p~9[| , except that they 
are built from scalar particles (spin-0) instead of electromagnetic fields, i.e. spin-1 bosons. If scalar fields exist in 
nature, such localized soliton-type configurations kept together by their self-generated gravitational field can form 
within Einstein's general relativity. 

In the case of complex scalar fields, an absolutely stable branch of such non-topological solitons with conserved 



particle number exists. In the spherically symmetric case, we have shown via catastrophe theory [|20 21 that these 
boson stars have a stable branch with a wide range of masses and radii. 

Kaup's first investigation ^] of the spherically symmetric boson star (BS) was based on massive scalar particles. 
Lateron, a nonlinear [/(|$| 2 ) potential was introduced by Mielke and Scherzer |23|], where also solutions with nodes, 
i.e. "principal quantum number" n > 1, were found. In building macroscopic boson stars, a nonlinear Higgs type 
self- interaction potential t/d^l 2 ) was later considered j24J] as an additional repulsive interaction. Thereby the Kaup 
limit for boson stars can even exceed the limiting mass of 3.23 Mq for neutron stars |25|] . 

Three surveys [^6|-p8| summarize the present status of the non-rotating case, a more recent survey including the 
rotating BS can be found in |29]. 

Recently, we construct f30| , |31|| for the first time the corresponding localized differentially rotating configurations via 
numerical integration of the coupled Einstein-Klein-Gordon equations. Due to gravito-magnetic effect, the ratio of 
conserved angular momentum and particle number turns out to be an integer a, the azimuthal quantum number of 
our soliton-type stars. The resulting axisymmetric metric, the energy density and the Tolman mass are completely 
regular. 

E. Are fundamental scalar fields part of nature? 

The physical nature of the spin-O-particles out of which the boson star (BS) consists is still an open issue. Until 
now, no fundamental elementary scalar particle has been found in accelerator experiments which could serve as the 
main constituent of the boson star. In the electroweak theory of Glashow, Weinberg, and Salam, a Higgs boson dublctt 
(<I> + , $°) and its anti-dublett (<E>~, $°) are necessary ingredients in order to generate masses for the and Z° gauge 
vector bosons [^2| . After symmetry breaking, only one scalar particle, the Higgs particle h :— (<E>° + $°)/v / 2, remains 
free and occurs as a state in a constant scalar field background. Nowadays, calculations of the two-loop electroweak 
effects enhanced by powers of the mass of the rather heavy top quark |53| has lead to an indirect determination of 
the Higgs mass, cf. @. For M t = 173.8 ± 5 GeV/c 2 , one finds m h = 104^ GeV/c 2 . So far, the experimental 
constraints are weak; even for the unrealistic case of a Higgs mass up to 1000 GeV/c 2 , the discrepancies for, e.g., the 
mass of the W boson would be small. Above 1.2 TeV/c 2 , however, the self-interaction U($) of the Higgs field is so 
large that the perturbative approach of the standard model becomes unreliable. Therefore a conformal extension of 
the standard model with gravity included could be necessary, see |3j|. Fermilab's upgraded tevatron has a mass 
reach of 135 < mh < 186 GeV/c 2 , while the high-energy experiments at the LHC at Cern will ultimately reveal if 
these Higgs particles really exist in nature. 

As free particles, the Higgs boson is unstable with respect to the decays h — > W + + W~ and h — ► Z° + Z°. In 
an hypothetical compact object like the BS, these decay channels are expected to be in equilibrium with the inverse 
pro cess Z° + Z° — > h, for instance. Presumably, this is in full analogy with the neutron star [ [37| or quark star 
|38]-[40|], where one finds an equilibrium of j3- and inverse /3-decay of the neutrons or quarks and thus stability of 



the macroscopic star with respect to radioactive decay. Such Higgs sector nontopological solitons |41| may also be 
candidates for cold dark matter. Nishimura and Yamaguchi [^2| constructed a neutron star using an equation of state 
of an isotropic fluid built from Higgs bosons. 

Nowadays there are many attempts of unifying the standard model with gravity on the quantum level, like string 



theory p3| . Commonly, the four-dimensional effective models make the prediction 44 that the tensor field g^ v of 
gravity is accompanied by one or several scalar fields. 

In string effective supergravities Eq] , the mass of the dilaton ip can be related to the supersymmetry breaking scale 
tosusy by m v ~ 10~ 3 (msusY/ TeV/c 2 ) 2 eV/c 2 with interesting astrophysical implications |Q, but this is not the 
only possibility. 
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The other scalar field of the effective string Lagrangian is the 'universal' invisible axion a, a pseudo-scalar potential 



for the Kalb-Ramond three form H := e v *da. (There are some speculations 47|4q| to identify it with the the axial 
part of a possible torsion of spacetime). From the Hubble scale H cq ~ 10 -27 eV/c 2 of matter-radiation equilibrium 
and the temperature T m ~ 100 MeV of mass generation at the epoch of chiral symmetry breaking, one can derive 
HQ the condition m a > (T m /eV) 2 H cq . This allows a very light axion mass m a = 7.4 x (10 7 GeV//cr) eV/c 2 
> 10 -11 eV/c 2 with decay constant f a close to the inverse Planck time, thus a prime candidate for dark matter. (This 
should not be confused with the Goldstone boson a of the Peccei-Quinn symmetry of standard QCD, for which a 
recent experiment | p2| has excluded the range of m a ~ 10~ 6 eV/c 2 . From cooling neutron stars, there can be inferred 
p3| an upper limit of m a < 0.06 - 0.3 eV/c 2 , depending on the equation of state of the nucleon fluid. 

II. BOSON STARS 

In a 1968 perspective paper, Kaup |2^] has studied for the first time the full generally relativistic coupling of linear 
Klein-Gordon fields to gravity in a localized configuration. This 'Klein-Gordon geon' is nowadays christened mini- 
boson star and can be regarded as a macroscopic quantum state. It was already realized that no Schwarzschild type 
event horizon occurs in such numerical solutions. Moreover, the problem of the stability of the resulting scalar 'geons' 
with respect to radial perturbations is treated. It is shown that such objects are, below a well-defined critical mass, 
resistant to gravitational collapse. These considerations are on a semiclassical level, since the Klein-Gordon field is 
treated as a classical field. 

The Lagrangian density of gravitationally coupled complex scalar field $ reads 

cbs = 4?^ + Kfertwca,*) - u(\ $ i 2 )]} , (2.i) 

where K = 8ttG is the gravitational constant. Using the principle of variation, one finds the coupled Einstein-Klein 
Gordon equations 

G„„ := - ^g^R = -/sT MV ($) , (2.2) 
a + 7^V=0, (2-3) 



where 



d I $ 



2 

is the stress-energy tensor and □ := ( g \ ) ( yj\ g \g^ v d v \ the generally covariant d'Alembertian 



A. Spherically symmetric solutions 



The stationarity ansatz 



$(r,t) = P(r)e" 



describes a spherically symmetric bound state of the scalar field with frequency u>. 
In the case of spherical symmetry and isotropic cordinates, the line-element reads 



ds 2 



3^dt 2 



dr 2 + r 2 (d6 2 + sin 2 9dy 2 \ 



(2.5) 



(2.6) 



in which the functions v = v(r) and A = A(r) depend on the radial coordinate r 



a/Tc 2 



+ y z + 



In a nut-shell, a boson star is a stationary solution of a (non-linear) Klein-Gordon equation in its own gravitational 
field; cf. |54|,|55|]. As in the case of a prescribed Schwarzschild background pq|, the spacetime curvature affects the 
resulting radial Schrddinger equation 



[d„2 - V eS (r) +uj 2 -m 2 ]P = 



(2.7) 
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for the radial function P(r) essentially via an effective gravitational potential V e s(r), when written in terms of the 
tortoise coordinate dr* :— e^ x ~ v ^ 2 dr . Then it can be easily realized that localized solutions fall off asymptotically 
as P(r) r-j (1/r) exp [—y/rn? — u> 2 r) in a Schwarzschild-type asymptotic background. 

The energy-momentum tensor becomes diagonal, i.e. T^ u (<&) — diag (p, — p r , —p_\_, —px.) with 

p=\(u?P 2 e- v + P' 2 e- X + U) , 

Pr = P-U , 

rl = Pr - P' 2 er x . (2.8) 

This form is familiar from fluids, except that the radial and tangential pressure generated by the scalar field are in 
general different, i.e. p r 7^ p±, due to the different sign of (P') 2 in these expressions. 

In general, the resulting system of three coupled nonlinear equations for the radial parts of the scalar and the 
(strong) gravitational tensor field has to be solved numerically. (Exact solution of massless scalar fields |37j or of the 
coupled Maxwell-Einstein-Klein-Gordon equation j5f| tend to be plagued with singularities.) 

In order to specify the starting values for the ensuing numerical analysis, asymptotic solutions at the origin and at 
spatial infinity are instrumental. The resulting configuration turns out to be completely regular and does not exhibit 
an apparent event horizon, cf. j3f). 

The stress-energy tensor of a BS, unlike a classical fluid, is in general anisotropic ]22| . In contrast to neutron 
stars 59; 6^], where the ideal fluid approximation demands an isotropic symmetry for the pressure, for spherically 
symmetric boson stars there are different stresses p r and p± in radial or tangential directions, respectively. The 
fractional anisotropy aj := (p r — p±_)/p r = P' 2 e~ x /(p — U) depends essentially on the self-interaction; cf. Rcf. J6l| , |62[ . 
So, the perfect fluid approximation is inadequate for boson stars. 

There exists a decisive difference between self-gravitating objects made of fermions or bosons: For a many fermion 
system the Pauli exclusion principle forces the typical fermion into a state with very high quantum number, whereas 
many bosons can coexist all in the same ground state (Bose-Einstein condensation). This also reflects itself in the 
critical particle number N := j°dr of stable configurations: 

• iVcrit — (Afpi/m) 3 for fermions 

• N CI - lt ~ (Afpi/m) 2 for massive bosons without self-interaction. 

Cold mixed boson-fermion stars have been studied by Henrique et al. p3[ and Jetzer [E6[. 



B. Critical masses of boson stars 



Since boson stars are macroscopic quantum states, below a certain critical mass M cr it they are prevented from 
complete gravitational collapse by the Heisenberg uncertainty principle ArAp ~ nh, cf. Ref. fl64f . This provides us 
also with crude mass estimates: For a boson to be confined within the star of effective radius i? e ff := (V-^0 j°rdr, 
the Compton wavelength of the collective boson has to satisfy A$ = (2n%/mc) < 2R e g. On the other hand, the 
star's radius should be of the order of the last stable Kepler orbit 3i?s around a black hole of Schwarzschild radius 
i?s := 2GM/c 2 in order to avoid an instability against complete gravitational collapse. 

For a mini-boson star, i.e. a massive boson model with merely the mass term [/(|$|) = m 2 |$| 2 as self-interaction, 
of an effective radius R c g = (ir/2) 2 Rs close to the last stable Kepler orbit of a black hole, one obtains the estimate 

M crit = (2/ 7 r)M| 1 /m > 0.633 M^/m - M Kaup , (2.9) 

cf. Ref. [^6|,^9). This provides us with a rath er g ood upper bound on the so-called Kaup limit. The correct value in 
the second expression was found numerically ]22fl as a limit of the maximal or critical mass of a stable mini-BS. Here 
Mpi := \JhcjG is the Planck mass and m the mass of a bosonic particle. 

For the likely mass Wh = 100 GeV/c 2 of the Higgs particle, e.g., one can estimate the total mass of this mini-boson 
star to be M ~ 10 10 kg and its radius R c s — 10~ 18 m yielding a density 10 45 times that of a neutron star. 

A boson star is an extremely dense object, since non-interacting scalar matter is very "soft". However, these 
properties are changed considerably by considering a repulsive self-interaction 

U(m = ™ 2 m 2 (l + |a(|<1>| 2 )^ = m 2 on |$ rcn | 2 . (2.10) 
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where A(|<I>| 2 ) denotes an arbitrary nonlinear self-interaction. The choice A(|$| 2 ) = (4A/m 2 )|<i>| 2 would lead us back 
to the quartic self-interaction of Colpi et al. p3| . If we adopt the value |<f>o| — Mpi/v'lSTr inside the boson star, one 
finds for the energy density 



p~m 2 |$ | 2 (l + A/8) 



(2-11) 



where A := A(Mp 1 /167r) is a dimensionless coupling constant such that we would recover A := (XMp^/Airm 2 ) for the 
quartic interaction. The self-interaction becomes dominating for A > 8 or A > 327r(m/Afpi) 2 . Thus, even a rather 
tiny coupling constant A could have drastic effects on a BS. 

Formally, this corresponds to a star formed from non-interacting bosons $ re n = <5(l+A/8) with a lower renormalized 
mass m — > m ren := mj \J\ + A/8 but larger Compton wave length Aj£ ren ) and, consequently, a larger radius R e s- (A 
reverse rescaling of the mass, as was presumed in a recent pre print p5| , leads to a smaller Compton wave length and 
other inconsistencies.) Consequently, we can apply again ( |2.9| ) for the Kaup limit and find that the maximal mass of 
a stable BS scales approximately as 



M r 



(2/ 7 r)M| 1 /m r 



M 2 

A/8^£l 



M 2 
A ^P1 



for A 



(2-12) 



For the quartic self- interaction, this accounts rather well for the numerical results of Colpi et al. P4| . Our formula 
not only reprod uces their asymptotic mass formula (11) for A — > oo, but, by construction, interpolates as well with 
the Kaup limit (Q for A = 0. 



Compact 
Object 


Critical mass 

M cri t 


Particle Number 

A crit 


Fermion Star: 
Mini-BS: 
Boson Star: 
Soliton Star: ||,|7) 


Mch := M^/m 2 
M K au P = 0.633 M 2 Jm 

10- 2 {M^/m^D 


~ (M P1 /m) 3 
0.653 (Mpi/m) 2 
~ (Mpi/m) 3 
2 x 10- 3 (M^/m 2 $g) 



1.5(GeV/mc 2 ) 2 M Q , where Mq denotes the solar mass. In astro- 
OmVXMlJm 2 = O.lVA (GeV/mc 2 ) 2 M Q which for A = 1 and 



The Chandrasekhar limit is Mch := Mp^/m 2 
physical terms, the maximal BS mass is M cr jt 

proton mass m ~ 1 GeV/c 2 is in the interesting mass range ~ 0.1 M Q of MACHO's. 

In a scale-invariant theory built from nonlinearly coupled dilatons tp, there arise a conserved dilaton charge via 
Noether's theorem from Weyl rescaling and thus will e nsure the stability of the configuration. For a dilaton star with 
quadratic self-interaction |67t , the same formula (2.12) applies, but the coupling constant A := (\M/4ttuj) 2 will be w 
dependent. For a very light dilaton <p of mass mdii = lO^ 11 eV/c 2 , resembling a misaligned 'universal' axion at its 
lower mass bound, Gradwohl and Kalbermann fl68|] found 
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VA 



Mr. 



rit — 40\/A km , 



(2-13) 



where A := \(M/ui) 2 is the rescaled coupling constant of the ip 4 interaction. 

To repeat, in building macroscopic boson a Higgs-type self-interaction f (|$|) is crucial for accommodating a 
repulsive force besides gravity. This repulsion between the constituents is instrumental to blow up the boson star so 
that much more particles will have room in the confined region. Thus the maximal mass of a BS can reach or even 
extend the limiting mass of 3.23 Mq for neutron stars ]2^ , |69| , |70| with realistic equations of state p = p(p) for which 
the (phase) velocity of sound is v s — ^Jdp/dp < c. However, this fact depends on the strength of the self-interaction. 

Therefore, if scalar fields would exist in nature, such compact objects could even question the observational AGN 
black hole paradigm in astrophysics. 



III. EXCITED BOSON STARS 



A. Gravitational atoms as boson stars 



Ruffini and Bonazzola [pi 



71 1 used the formalism of second quantization for the complex Klein-Gordon field and 
noticed an important feature: If all scalar particles are within the same ground state |$) = (N,n,l,a) — (TV, 0,0,0), 
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which is possible because of Bose-Einstein statistics, then the semi-classical Klein-Gordon equation of Kaup can 
be recovered in the Hartree-Fock approximation for the second quantized two-body problem. In contrast to the 
Newtonian approximation, the full relativistic treatment avoids an unlimited increase of the particle number N and 
negative energies, but induces critical masses and particle numbers with a global maximum. 

Due to this Hartree-Fock approximation and while als o neglecting effects of the quantized gravitational field, the 



same coupled Einstein-Klein-Gordon equations ( p . 2|j2 .3|) apply. Therefore, a boson star is also called a gravita- 
tional atom [Jr2[ . Since a free Klein-Gordon equation for a complex scalar field is a relativistic generalization of the 
Schrddinger equation, we consider for the ground state a generalization of the wave function 

\N, n,l,a): 3 = K(r) Y^O, ^e^'^ 

' R n a {r)P?{cos6)e laip e- l{E ~ /h)t (3.1) 



-in 

of the hydrogen atom. Here R™{r) is the radial distribution, Y l a (9, tp) the spherical harmonics, P ; a (cos#) are the 
normalized Legendre polynomials, and \a\ < I are the quantum numbers of azimuthal and angular momentum. 

Due to their inherent 'gravitational confinement' gravitational atoms represent coherent quantum states, which 
nevertheless can have macroscopic size and large masses. The gravitational field is self-generated via the energy- 
momentum tensor, but remains completely classical, whereas the complex scalar fields are treated to some extent as 
Schrddinger wave functions, which in quantum field theory are referred to as semi-classical. 

Moreover, Feinblum and McKinley |73| found eigensolutions with nodes corresponding to the principal quantum 
number n of the H-atom. Motivated by Heisenberg's non-linear spinor equation |74| additional self-interacting terms 
describing the interaction between the bosonic particles in a "geon" type configuration were first considered by Mielke 
and Scherzer pj^ j, where also solutions with nodes, i.e. "principal quantum number" n > 1 and non- vanishing angular 

momentum I ^ for a t'Hooft type monopole j7fj] ansatz <& I=a ~ R{r) Pj a ^ (cos 9) were found. These highly interesting 
instances of a possible fine structure in the energy levels of gravitational atoms poses the question if quantum geons 
jrF|p"9[ are capable of internal excitations? Recently, without reference to these earlier works, such "exited boson stars" 
were recovered |7(],|77j and their stability properties corroborated in some numerical details. Moreover, Rosen |7§| ] 
reviewed his old idea of an elementary particle built out of scalar fields within the framework of the Klein-Gordon 
geon (or the mini-boson star, as they are christened today). 

Several surveys JUJIEIzEeI summarize the present status of the non-rotating case. 



B. Rotating boson stars 



In the framework of Newtonian theory, boson stars with axisymmetry have been constructed by several groups. 
Static axisymmetric boson stars, in the Newtonian limit [8C| l and in GR |8]J| , show that one can distinguish two 
classes of boson stars by their parity transformation at the equator. In both approaches only the negative parity 
solutions reveal axisymmetry, while those with positive parity merely converged to solutions with spherical symmetry. 
The metric potentials and the components of the energy- momentum tensor are equatorially symmetric despite of the 
antisymmetry of the scalar field. In the Newtonian description, Silveira and de Sousa |82| followed the approach of 
Ref. (7^] and constructed solutions which have no equatorial symmetry at all. Hence, in GR, we have to separate 
solutions with and without equatorial symmetry. 

Kobayashi et al. H?) tried to find slowly rotating states (near the spherically symmetric ones) of general relativistic 
boson stars, but they failed. The reason for that is a quantization of the total angular momentum [ |S4j 

J = / T 3 ° ^/\g\d 3 x = aN a = 0, ±1, ±2, • • • (3.2) 

of boson stars which is proportional to the particle number N and vanishes if a = 0. This relation between angular 
momentum and particle number was first derived by Mielke and Schunck |3l|] . 

In recent papers ]30| , |85t , we proved numerically that rapidly rotating boson stars with a / exist in general 
relativity. Because of the finite velocity of light and the infinite range of the scalar matter within the boson star, 
our localized configuration can rotate only differentially, but not uniformly. This new axisymmetric solution of the 
coupled Einstein-Klein-Gordon equations represent the field-theoretical pendant of rotating neutron stars which have 
been studied numerically for various equations of state and different approximation schemes [[39 
for differentially rotating superfluids |8q | as a model for (millisecond) pulsars. 

On the basis of Ref. HI, it has erroneously been claimed (8lj that "rapidly rotating boson stars cannot exist". 
However, more recently |89[] the same Japanese group (as well as Q) followed exactly our Ansatz and could verify 



70 



as well as 
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all our earlier results pd,Bl[, albeit of some extension to stronger gravitational fields, due to better computational 
facilities. Due to the anisotropy of the stress-energy tensor, our configuration is differentially rotating, see i fu) for 
more details. 

Moreover, the energy density of our rotating boson star is concentrated in an effective mass torus |p5| . Thus this 
first nonsingular model of a rotating body in GR realizes to some extent the suggestion of Newman et al. J9l|] to fill 
in the Kerr metric, in view of its ring singularity, with a toroidal rather than a spherical source. Toroidal structure 
occurs also in relativistic star systems with an accretion disk |9S|. 

Since rotating BS have a toroidal structure, there seems to exist the more speculative possibility of knotted vortex 
like excitations, cf. Ref. |B3[. For an 0(3) Skyrme model, their existence has recently been demonstrated numerically 




C. Formation of (primordial) boson stars 



The possible abundance of solitonic stars with astrophysical mass but microscopic size could have interesting 
implications for galaxy formation, the microwave background, and formation of protostars. The formation of non- 
gravitating non-topological solutions was already studied by Frieman et al. |9(J.) In comparison with primordial black 
holes, it is therefore an important question if boson stars can actually form from a primordial bosonic "cloud" |97|. 

Collisionless star systems are known to settle to a centrally denser system by sending some of their members to 
larger radius. Likewise, a bosonic cloud will settle to a unique boson star by ejecting part of the scalar matter. 
Since there is no viscous term in the KG equation ( |2.3| ), the 'radiation of the scalar field is the only dissipationless 
relaxation process called gravitational cooling. Seidel and Suen |9sU99f| demonstrated this numerically by starting 
with a spherically symmetric configuration with Mi n i t i a i > Mxaup, i-e- which is more massive then the Kaup limit. 
Actually such oscillat ing an d pul sating branches have been predicted earlier in the stability analysis of Kusmartsev, 
Mielke, and Schunck [100,21,101 by usi ng c atastrophe theory. Oscillating soliton stars were constructed by using real 
scalar fields which are periodic in time 102]. Without spherical symmetry, i.e. for $ ~ R a (r)Yi a (6 , ip), the emission 
of gravitational waves would also be necessary. 

For the formation of primordial BSs, an important issue is the breaking of unified gauge (super-)symmetry at 
high temperature in order to yield a scalar-antiscalar asymmetry e$ = N^/N^, as in the case of baryon-antibaryon 
asymmetry, where €b = Ns/N-y ~ 10 -10 . Here, we recall the situation of collapsing homogeneous mini-BS clouds in 



the early Universe, cf. [103,281. Because the Jeans scale at decoupling time is greater than the horizon scale, a bosonic 



mass of Afpi/m 2 immediately collapses and since this is a factor Mp\/m higher than the maximal mass allowed within 
the mini-BS model, only black holes can form. For an asymmetry factor of order es ~ ro/Mpi, however, the total 
mass remaining within the horizon is Al^/m, hence, BSs could form, avoiding the final state of a black hole. 

For a real (pseudo-) scalar field, like the axion a of the broken Peccei-Quinn symmetry J5lj in QCD, the outcome 
is quite different. The axion has the tendency to form compact objects (oscillatons) in a short time scale. Due to its 
intrinsic oscillations it would be unstable, contrary to a BS. Since the field disperses to infinity, finite non-singular self- 
gravitating solitonic objects cannot be formed with a massless Klein-Gordon field [104,105|. In Ref. [ 106 1 a different 
mechanism for forming axion miniclusters and starlike configurations was proposed. The self-coupling relaxion time 
[ 37| is compatible or larger as the age of the Universe. For fermionic soliton stars, there is a temperature dependence 
[ lOTf in the forming of cold configurations. 



D. Gravitational waves 



In the last stages of boson star formation, one expects that first a highly excited configuration forms in which the 
quantum numbers n, I and a of the gravitational atom, i.e. the number n — 1 of nodes, the angular momentum and 
the azimuthal angular dependence e lav are non-zero. 

In a simplified picture of BS formation, all initially high modes have eventually to decay into the ground state 
n = Z = <z = 0bya combined emission of scalar radiation and gravitational radiation. 

In a Newtonian approximation of Ferrell and Gleiser |72j , the energy released by scalar radiation from states with 
zero quadrupole moment can be estimated by 

£rad ~ (n - l)MlJm . (3.3) 

This is accompanied by a loss of boson particles with the rate AN ~ (n — l)(Mpi/m) 2 . For investigating the 
gravitational radiation of macroscopic boson stars with large self-interaction, a reduction of the differential equations 
can be taken into account [BO] . 
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The lowest BS mode which has quadrupole moment and therefore can radiate gravitational waves is the 3d state 
with n = 3 and I = 2. For AJ = 2 transitions, it will decay into the Is ground state with n — 1 and I = while 
preserving the particle number N. The radiated energy is quite large, i.e., E Ta ^ — 2.9 x 10 22 (GeV/mc 2 ) Ws. Thus 
the final phase of the BS formation would terminate in an outburst of gravitational radiation despite the smallncss 
of the object. 



E. Gravitational evolution 



There would occur an evolution of boson stars if the external gravitational constant n changes its value with time 
[ 108 - 1 12 1 . This can be outlined within the theory of Jordan-Brans-Dicke or a more general scalar tensor theory. 
The results show that the mass of the boson star decreases due to a space-depending gravitational constant, given 
through the Brans-Dicke scalar. The mass of a boson star with constant central density is influenced by a changing 
gravitational constant. Moreover, the possibi lity of a gravitational memory of boson stars or a formation effect upon 
their surrounding has been analyzed as well [111]. 



IV. ARE MACHOS AXION STARS? 



Direct observation of boson stars seems t o be impossible also in the far future. We propose here several effects 
which could possibly give indirect evidence [ 1 13 1 . In the asymptotic region, the rotation velocity of baryonic objects 
surrounding the boson star can reveal the star's mass. Assuming that the scalar matter of the BS interacts mainly 
gravitationally, we would have a 'transparent' BS detecting a gravitational redshift up to values of z = 0.68 observ- 
able by radiating matter moving in the strong gravitational po tenti al. For further investigations of rotation curves, 
cf. Ref. [114,115], and of boson stars as gravitational lenses, cf. [116]. 

Solutions with an infinite range can be found where the mass increases linearly |l5"|,[l6| . In the context of the dark 



matter hypothesis, it may be speculated if such boson halos as well as excited BS states 1 115 114] can be used to 
fit the observed rotation curves for dwarf and spiral galaxies [[l5 16 . Boson halos have a finite radius if a positive 
cosmological constant exists as most recent results from supernovae reveal ][|||]. 

Moreover, BSs could be the solution for the MACHO problem, as we are going to analyze in more detail. 

In effective string theories, the dilaton ip, another moduli field (3, and the 'universal' invisible axion a are predicted 
B| . This can be read off from the effective string Lagrangian 



R + g» v dn,<pd v <p - W*PB v p - -e 2 «^acU 



2A 



(4.1) 



This corresponds to Eq. (11) with 77 = 2 of Ref. [117] and allows to combine [ 1 18 1 the axion and the dilaton into a singl 



complex scalar field $ := a + ie v , the axidilaton. In the conformally related Einstein frame g M „ — > g M „ := e 
and for constant modulus j3, our results on BSs can easily by transferred to this axidilaton content of strings. 



9^ 



A. Mass range of axion stars 

As macroscopic quantum states, BSs arc quite generally prevented from complete gravitational collapse below a 



critical total mass M cr it which, typically, depends inversely on the particle mass, see Eq. (2.9). 

The numerical results are shown in Fig.[j]. The left figure exhibits the dependence of the mass M and the particle 
number TV (rest mass) on the central density po. Stable axionic BSs exist at central densities lower than the maximum 
mass. The critical values are: M cr it = 0.846 M Q , miVcrit = 0.873 M© and p c = 9.1 x p nuc \, where p nuc \ = 2.8 x 10 17 
kg/m 3 is the average density of nuclei. Since non-interacting bosons are very "soft" , BSs are extremely dense objects 
with a critical density higher than for neutron or strange stars p0| . The figure on the right hand side gives the mass 
depending on the radius (measured in km). For the mass-radius diagram, we have chosen as radius 99.9% of the 
total mass. This ensures that the exponentially decreasing 'atmosphere' of the BS has almost no influence on the 
asymptotic Schwarzschild spacetimc. 

The stable BSs or axion stars (ASs) have radii larger than the minimum at 20.5 km and a mass of 0.846 M@. In 
order to derive these values, we have assumed that the mass of the scalar field is 10 -10 eV close to the lower bound 
of axions, leading to an asymmetry factor of eg ~ 10 -38 . and that no self- interaction exists. We stress that the 
total mass of these relativistic ASs is just in the observed range of 0.3 to 0.8 M© for MACHOs. One could also turn 
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this argument around: By identfying the MACHOs with known gravitational mass of about 0.5 M Q with ASs, we 
are essentially 'weighing", via MKaup/A^rit — m i the axion mass to m a ~ 10 -10 eV/c 2 . It is gratifying to note that 
such a low value is perfectly compatible with the constraints on the mass range of the Kalb-Ramond axion seeding 
the large-scale CMB anisotropy, cf. the recent results of Gasperini and Veneziano (5(]^j| within low-energy string 
cosmology. 



For the other option of dilaton ip being stabilized [ 11£ | through the axion, a much smaller dilaton mass of 



m 



10~ 6 m a could be generated non-perturbatively, such that the dilaton behaves very similar to misalignment produced 
Peccei-Quinn axion a. Our conclusion also with respect to the mass range of an axidilaton star will not changed 
much, if we use the full Brans-Dicke type interaction |68fl for the combined axidilatons. 



Thus, for cosmologically relevant (invisible) axions as cosmological dark matter also an AS [97,113 106 1 with a 
rather large mass of would be possible and stable. 

Therefore, if such-string inspired scalar fields would exist in Nature, axions could not only solve the non-baryonic 
dark matter problem ||, but their gravitationally confined mini-clusters, the axion stars, would also represent the 
observed MACHOs in our Galaxy. 

V. OUTLOOK: GRAVITATIONALLY CONFINED HAWKING RADIATION? 

Commonly for the Bekenstein-Hawking radiation the spacctimc geometry is treated as a given fixed background, 
e.g. the Schwarzschild solution. However, due to the universality of gravitational interaction, the evaporating quantum 
field, say a scalar field may have a "back-reaction" upon the spacetime geometry via the semi-classical Einstein 
equation 

= — k(0|T m1/ ($)|0) . (5.1) 



For instance, a 'bouncing shell' model [120| with retarded time u leads to (0|T Mtl |0) — > k 2 /487t, the standard Hawking 
result. The situation becomes, however, much more complicated by the fact that the vacuum expectation value 
(0 1 0) of the energy-momentum tensor T^ V1 for instance defined by the point-splitting prescription, is not unique. 
One ambiguity in (0|T /JJ/ ($)|0) is of the type rn^G^, i.e. linear in the curvature, and can be readily absorbed in a 
redefinition of the 'bare' gravitational constant k. However, the next order corrections are quadratic in the curvature 
and therefore of the same one-loop order arising from the notorious nonrenormalizability of perturbative quantum 



gravity, cf. Ref. [ 1 2 1 1 , p. 90. To some extent, the finite part of such higher order curvature counterterms in the 
Lagrangian can be simulated by a self- interaction potential U($), cf. ]122|] . 

Already on the semiclassical level one could ask the question what happens to the (massive) particles associated 
with the second quantized field $ in a patch of some strong gravitational background field? Could the particles created 
by the Unruh effect instead of evaporating to infinity rather form a bound state within their self- consistently generated 
gravitational field? Moreover, could it be possible that the full back-reaction on the geometry is strong enough lead to 
a curved spacetime without horizon and singularities, similarly as in some exact solvable (2+l)-dimensional models? 
Actually some aspects of this issue were already answered by Ruffini and Bonazzola |6l|] for a spherically symmetric 



self-gravitating configuration of N particles in a Hartree-Fock approximation. Thus the back-reaction (5.1) may lead 
us back exactly to some stable branch of boson stars where the particles are treated on the first quantization level. 
These type of stars have an exponentially decreasing energy density of the scalar field, an 'exosphere ' of particles in 
the stable state of equilibrium of particle creation and annihilation. Moreover, for these type of compact objects wit h 
an effective radius close to the last stable Kepler orbit an event horizon is suppressed due to the back-reaction ([} . if) . 

Below the Kaup limit, we have seen that such macroscopic quantum states are absolutely stable, at higher central 
densities the configuration becomes more and more unstable, and undergoes complete gravitational collapse. 

So could it be that the picture of an evaporating black hole is just a first order semi-classical approximation; rather, 
below some mass limit, we may end up in a self-consistent state of a boson or fermion star with a gravitationally 
confined Hawking radiation, a quantum geon? 
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FIG. 1. Left: Mass M and particle number N (or rest mass mN at infinity) of a BS depending on the central density p. 
Right: Mass-radius dependence of an axionic BS. 
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